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The anomalous exponent, r/p, for the decay of the reunion probability of p vicious 
walkers, each of length A^, in d (= 2 — e) dimensions, is shown to come from the 
multiplicative renormalization constant of a p directed polymer partition function. 
Using renormalization group(RG) we evaluate r]p to O(e^). The survival probability 
exponent is rjp/2. For p = 2, our RG is exact and rjp stops at 0(e). For d = 2, the 
log corrections are also determined. The number of walkers that are sure to reunite 
is 2 and has no e expansion. 
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A question of perennial interest is this: If p vicious (i.e. mutually avoiding) random 
walkers in d dimensions start from a point in space at time z = 0, how does the probability 
of reunion at time z = N decay with A^? A complementary question is the decay of the 
survival probability. Generally, for large A^, a power law form is expected signifying certain 
universality in the behaviour which we explore in this paper. 

One of the main motivations to study vicious walkers is to understand and predict the 
nature of the phase transitions for interface wetting phenomena, melting of the commensu- 
rate phase, commensurate - incommensurate (CI) transitions etc, all of which involve walls 
or string like objects These walls, which are statistically parallel to the z axis, with 

fluctuations in the transverse direction may reunite if defects are present in the system. The 
loops formed thereby are often relevant for criticality . They are also important in 

the context of abelian sand pile model of self organized criticality. These directed strings 
can be embedded in arbitrary d + 1 dimensions with d as the transverse spatial dimension. 
These are called directed polymers (DP). Note that the preferred direction of the DP plays 
the role of time for the random walker. 

In spite of its enormous applications, exact results (e.g., the decay exponents) available 
for reunion of a system of vicious random walkers are restricted to d = 1 dimension only, 
with an approach for higher dimensions remaining as a challenge. It is clear that the 
constraints for the walkers to be reunited at one or both the ends and their repulsion do 
change the exponents drastically. Our focus is on these exponents for d > 1. 

In this paper, we adopt the DP description. The viciousness becomes an equal time 
mutual repulsion for the DPs. The reunion and the survival probablities then come from 
appropriate partition functions for these interacting DPs. Since we are interested in the 
asymptotic behaviour, i.e., in the large chain length limit, we adopt a continuum model. 
The merit of the path integral technique easily enables us to formulate the model and obtain 
the exponents of the chain length for the partition functions in arbitrary d + 1 dimensions 
through a renormalization group (RG) approach. In this formulation, the dimensionless 
hamiltonian for p DPs in d + 1 dimensions, each of length A^, is given by (Note: z is the 
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preferred direction) |]8|J^,p^ 

p N N 

= 2 E / i-? + ^0 E / dz5\v,,iz)) (1) 

i=l 

where rj = dri/dz, and rij{z) = ri{z) — rj{z), ri{z) being the d dimensional coordiante of a 
point at z on the contour of the ith chain. The first term is the elastic energy, taking care 
of the connectivity of the chains. The viciousness of the walkers is simulated in the second 
term by a 5-function interaction at the same z coordinate with fo > 0. Dimensional analysis 
shows that Vq is dimensionless and hence marginal at d = 2 (upper critical dimension). The 
fact that d = 2 is special will be reflected in the later discussion. Our method can be 



extended to many body interactions to be discussed elsewhere [O. 

The reunion probability for p chains follows from the constrained partition function Z^^p 
for DPs with the restriction that the chains are all tied together at both the ends at origin. 
A computationally easier quantity is the no-reunion or survival probability. This comes from 
the partition function, Zs^p, with the restriction that the walkers are all tied together at the 
origin intially, i.e., at 2; = but are free at the other end. These and the exponents are 
defined generically below {g = R or S) as 

Zg^p = jvvFg e-"^ ~ N-^^-". (2) 

Here, Vr stands for the sum over all paths, and Fg denotes the end point constraint, im- 
plemented by a 5'^(rj(0)) for one end of each chain {g = R, S) and a 6'^{ri{N) — r) for the 
other end of each chain for reunion at r. One can consider a more general reunion problem 
where the walkers can meet anywhere in space at 2; = A^. This, 2i?,p, requires an extra 
integration of the partition function of the Eq. |] type over the end point coordinate, and 
^R,p ~ N'^'^S'P. We have indeed checked explicitly for a few cases that, the anomalous part 
of ipji^p and "^R^p are the same. Incidentally, these exponents have recently been used to 
discuss a novel crossover in the Ising model 0. 

Once the exponents are known, one can address the question of the critical number 
of vicious walkers that are sure to reunite in d dimensions. Conventional analysis [see. 



e.g., mi]] shows that this number is determined by the convergence criterion of the conditional 
probabiUty Z^p/Zs^p ~ A^"-^, [x = ^_R,p — ^5,p], i-e., whether / dNN~^ is convergent or not. 
The critical number is then obtained from ^ R^p — ips,p = 1- We call this number pc- 

Our approach is to start with a perturbation expansion (digrammatically) in the coupling 
constant, and use the renormalization group (RG) approach to go beyond the validity of 
the perturbation series [^|lOl. The series stumbles on divergences that require standard 
dimensional regularization to identify the poles. The removal of the divergences demands a 
renormalization of the interaction parameter and an overall multiplicative renormalization 
constant [|l^]. This constant gives rise to an anomalous exponent of the length scale which 
is unexpected from a conventional dimensional analysis. Such constants were not needed for 
virial coefficient or other thermodynamic properties of DPs and hence no anomalous 



dimension ever appeared there. 



We, in this paper, first study the two chain reunion problem through an exact RG [jlS 
It is then extended to the p chain case. We also study in a similar way the no-reunion case, 
without giving much details. The results then produce a surprising "super" universality for 
Pc- Our procedure is general enough to answer many other questions than studied here. 
These will be discussed elsewhere. 

Had the walkers been noninteracting {vq = 0), the exponents would follow trivially from 
the normalized distribution G{r \ z) = (27rz)^^/2 exp(— r2/2z), for one walk of length z and 
the end to end distance r. The "gaussian" exponents are 

'4's,p = 0, V'R.p = Pc?/2, and ^ij,p = (p - l)d/2. (3) 

These gaussian exponents differ from the exact d = 1 results for interacting walkers 

i^s,p = ^ , ^R,p = y, and "^R^p = ^ . (4) 

The difference between any of these and the corresponding gaussian exponent is the anoma- 
lous dimension to be denoted by rjs^p (for survival) and rjp (for reunion). A consequence of 
Eqs. ^ and ^ is that, for d = 1, pc = 2 for vicious walkers in contrast to 3 (1 + 2/c/ for 
general d) for noninteracting walkers. 



One might expect pc = 2 for c? < 2 on a simple argument. For two vicious walkers, 
one can consider the relative walker who starts from one of the nearest neighbours of the 
origin (NNO) but avoids the origin. Since for = 2 a random walk is recurrent, one might, 
ignoring the origin avoidance, expect the walker eventually to come to one of the NNOs 
again, forcing a reunion. Our analysis shows that this is indeed the right answer. 

Two chain case: The simplest case that can be calculated exactly is the reunion problem 
of two walkers. The partition function which restricts the two walkers to be tied at the ends 
at spatial coordinate r = can be calculated with perturbation in the interaction parameter 
about the free walkers. The diagrams representing the terms of the perturbation series are, 
e.g., the first three diagrams of Fig la with the left most line omitted. The general term 
corresponding to DPs with arbitrary number of encounters can be calculated easily using 
the identity for m gaussian propagators 

G™(r I z) = (27rz)-("-i)'^/2m-'^/2G(r | z/m). (5) 

In terms of the dimensionless coupling constant Uq = VqL'^, {e = 2 — d and L an arbitrary 
length scale) the series for the partition function is 

oo 

{27rNfZn,2 = 1 + E(-«o/4vr)'^(47riVL-2)-/2^^^^(,/2). (6) 

n=l 

where Qn{x) = r"(x)/r(na;). This series, clearly showing divergences at e = 0, requires 
renormalization through absorption of these poles. We quote below, for convenience, the 
formula for renormalized u and the corresponding /3-function from Ref. [0 

Mo = u[l - u/{27re)]-\ and p{u) = ue [1 - u/{2ne)]. (7) 

with u = u* = 2'Ke as the fixed point where (3{u) = 0. This interaction renormalization is 
not sufficient to absorb even the divergence in the first order term in uq of Eq. ^. Let us 
define the renormalized partition function as 

Zr,pI = Rr,p{u) Zr^p withi?i?,p(M) = 1 + 51 (S) 

n)0 
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as the overall multiplicative renormalization constant {p = 2). The divergences can be 
absorbed completely provided bi = (yre)"^, and 62 = 3/(47r^e^). In fact, it is possible to sum 
the whole series to get Rr^2{u) = [1 — u/ {2-ne)Y'^ . 

An appeal to the renormalization group equation |]12| immediately shows that '4)r^2 is 



different from the naive engineering dimension (= d), Eq. ^by 772 = 7r,2(w*) obtained from 

d 

27i?,p(M) = /^(«)^ In RrAu) (9) 
with p = 2. We, therefore, have 

V2 = e, iIjr,2 = 2 and ^jj,2 = 2 - d/2. (10) 
in agreement with the d = 1 exact results of Eq. ^ |13 



p chains: Let us first consider the partition function Z^p for p chains tied together at 
the origin at both the ends. The terms in the perturbation expansion upto second order in 
Vq correspond to the diagrams shown in Fig. la. The zeroth order gaussian term (Fig. la.l) 
contributes {2nN)-P'^/'^. 

The second order diagrams are of several types: ladder type involving two chains (Fig. 
la. 3), three chains connected by interaction (Fig. la. 4), and two pairs connected indepen- 
dently (Fig. la. 5). The first and the last of these and the first order diagram (Fig.la.2) are 
already evaluated in Eq. since these are just two chain terms. The most crucial contri- 
bution is from Fig. la. 4 which connects three chains. The two chain terms are important for 
renormalization, but the O(e^) contribution to the exponents, ultimately, comes from this 
diagram. 

The contribution of Fig. la. 4 involves (2-7? A^)"^^"^)'^/^ for the noninteracting (p — 3) 
chains. To evaluate the connected piece, one has to use an identity similar to Eq. ^ for m 
gaussian propagators, and the markovian and normalization properties of G. Details will be 



given elsewhere ||rT|. The final expression is 

(27r7V)-P^/2iV^(4vr)-'^[6;3(e) 2^1 (e, e; 3e; 3/4) + 

(3/4)^ r(-e)r-i(l - e)^2(e) 3F2 (e, e, 1; 2e, 1 + e; 3/4) 



where e = e/2, 3F2 is the generahzed hypergeometric function, and ^„(x) is defined after Eq. 
^. Combining all the terms with appropriate symmetry factors (see Fig. Ic) and expanding 
each one around e = 0, we find that 

^""'^ -1-^0 (P^)(2 + e\nx) + -^(C~tD + tC\nx) + ... , (11) 



(27rAr)-P'i/2 27re ' An^e^ 

where C = Qip"^ -p+1), D = 3(f) ln(3/4), and x = AtxNL-^. 

Substituting tt, from Eq. 0, for uq (note that mq still pertains to two body 5 function 
interaction) in Z/jp, Eq. [Tl|, one sees that the series still requires an overall multiplicative 
renormalization constant Rr^p{u) = 1 + hiu + 62"^^ + ... (see Eq. |]). Minimal subtraction of 
the poles |12[ yields 61 = (D/(7re), and 62 = (C + e D)/{2-Kef. 



The anomalous exponent comes from ■jR^p{u) of Eq. ||. Using the /3-function, Eq. 0, we 

get 

27i?,p(^i) = + (262 - bi{2ne)-^ - bl)u + ..]. (12) 

Finiteness of 7_r,,p(m) in the limit e — is guaranteed by the exact cancellation of the 
0{e~^) term in 262 — 6i(27re)~^ — b1 leaving alone the 0(e~^) term. In fact, the finiteness 
criterion of 7_r,p(m) as e ^ dictates the cancellation of all but the 0(e^^) part in all higher 
orders of u. Therefore, to calculate the O(e^) term in 7('U*), no higher order term in the 
renormalization constant in Rr^p{u) is required. This O(e^) term can be traced back to the 
three chain connected diagram. 

We finally obtain from Eq. |12] the anomalous exponent for reunion as 

Vp = 7rA^*) = (^) e + 3 (?) ln(3/4) + 0(6^). (13) 

The survival case, can also be considered in an identical fashion. The diagrams are 
shown in Fig. 2b. The unexpected result we find is that, upto O(u^), the renormalization 
constant Rr^u) is the square of the renormalization constant needed for Zg^p- We have 



seen it exactly for p = 2. Though it seems rather obvious or intuitive |]T3, still a general 
proof is lacking. We conclude 



ips,p = Vs,p = Vp/'^- (14) 

This is a scaling relation not hitherto recognized at all. 

At the upper critical dimension d = 2, the /3-function has only the trivial fixed point. 
Its integration gives u{L) = uo[l + (Mo/27r) ln(L/Lo)]^"'^ where Uq = u{Lq). The renormal- 
ization group equation, furthermore, via the method of characteristics |T^, leads to the log 
correction for the reunion probability, as 

Zr,p ~ iV-(^-^) [uo ln{N/No)r'^^-'^/' (15) 

where A^^o = L^. Similarly, the survival probability, instead of remaining N independent as 
for the gaussian case, now has a slow deacy as [\n(N / Nq)]~p^p~^^^^ . 

To determine Pc, we now solve {p — l)d + r^p = 2 in an iterative way. Upto O(e^), Pc = 2, 
even after incluidng the log corrections at d = 2. Taking into account the exact result 
for d = 1, we predict that Pc = 2 for d < 2 and is independent of d. Of course, two walkers 
are not certain to reunite for d > 2. Still, the fact that the interaction can make the critical 
number "super" universal (at least weakly, i.e., for < 2), is surprising. 

In summary, we emphasize that for the problem of vicious walkers, that can be visualized 
via the path integrals as directed polymers, an e = 2 — d expansion is sucessful in obtaining 
the essential results. The key feature is the evaluation of the overall multiplicative renor- 
malization constants for the partition functions. The relevant exponents for p chains are 
computed to O(e^). Extensions to still higher orders in e seem not out of reach. Our analysis 
shows that for O(e^), three chains should be collectively aware of each other. Such collective 
features are expected to continue in higher orders also. For p = 2, we obtain exact 0(e) re- 
sults. A significant result is the scaling relation Eq. between the anomalous exponents for 
reunion and survival probabilities. The speciality of the dimension d = 2 shows up through 
the dependence of the probabilities on In A^ with p dependent powers. Furthermore, though 
3 (2) noninteracting walkers are sure to unite in one (two) dimensions, we find the number 
is 2 for vicious walkers, independent of d < 2. There seems to be a superuninversality about 



8 



this number. Numerical verification of these predictions, especially on fractals, would be 
welcome. 

A major part of this work was done at ICTP, Trieste, during our visit there under the 
Federation scheme. We thank ICTP for its warm hospitahty. 
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FIGURES 

FIG. 1. Diagrammatic expansion upto Vq for p > 3 chains tied at (a) both ends and (b) one 
end. The sohd hnes represent the chains and the dashed lines represent the interaction at same 
chain length z (the involved points on the chains have the same spatial coordinates.). Only a few 
chains are shown. Rules for evaluation: (i) a factor of G for each piece of the solid lines, (ii) —vq for 
each dotted line, (iii) an integration over the coordinates of the interaction points, and (iv) a time 
ordered integration of the z coordinates. See Ref. [p|,p!0|. (c) Symmetry factors. These numbers 
take into account the number of diagrams with identical contribution but generated either through 
a mere permutation of the chains or by reversal of time ordering. 
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